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In this paper, we investigate experimentally the non-equilibrium steady state of an active colloidal 
suspension under gravity field. The active particles are made of chemically powered colloids, showing 
self propulsion in the presence of an added fuel, here hydrogen peroxide. The active suspension 
is studied in a dedicated microfluidic device, made of permeable gel microstructures. Both the 
microdynamics of individual colloids and the global stationary state of the suspension under gravity 
- density profiles, number fluctuations - are measured with optical microscopy. This allows to 
connect the sedimentation length to the individual self-propelled dynamics, suggesting that in the 
present dilute regime the active colloids behave as 'hot' particles. Our work is a first step in the 
experimental exploration of the out-of-equilibrium properties of artificial active systems. 

PACS numbers: 



The collective behavior of 'active fluids', made of self- 
propelled entities, has raised considerable interest over 
the recent years in the context of non-equilibrium sta- 
tistical physics [THE]. Such systems are rather common 
in living systems, from swimming cells, bacteria colonies 
[9UT3], to flocks of birds or fishes [2]. These entities 
move actively by consuming energy and their behavior is 
thus intrinsically out of equilibrium. Building a general 
framework describing their collective properties remains 
accordingly a challenging task and led to a considerable 
amount of work towards this aim [IHH]. By contrast, 
much less work has been performed on the experimental 
side, and mainly on assemblies of living micro-organisms 
- which are naturally self-propelled - [TTHT4] , however at 
the expense of a lack of control and flexibility of the in- 
dividual particles and their interactions. There is there- 
fore a need for new experiments based on artificial model 
systems, involving suspensions of microscopic active par- 
ticles with controled propulsion and interaction mecha- 
nisms - here designated as 'active suspensions'. Several 
routes to design individual artificial microscopic swim- 
mers have been explored recently [T5HT7] , taking benefit 
of the recent progress made to shape and design colloidal 
particles at microscales. However going from the individ- 
ual to the many particles situation remains a challenging 
task and has not been achieved up to now with artificial 
motile particles. 

In this paper we explore experimentally the behavior of 
a dilute active suspension of articial swimmers under an 
external (gravity) field. This problem was discussed the- 
oretically in a recent contribution by Tailleur and Cates 
[8]. Here the active colloids are powered chemically, fol- 
lowing the route proposed by Howse et al. [15] . The asy- 
metric dismutation of hydrogen peroxide (H2O2) on the 
colloid itself is used as the driving power, on the basis of 
a self (diffusio-) phoretic motion, induced by the building 
up of an osmotic pressure gradient at the interface of the 
colloid [18-20 j. Furthermore a specifically designed mi- 




FIG. 1: (a) SEM picture of janus, platinum-latex, colloids 
(scale bar 500nm). (b) Sketch of the experimental setup: a 
circular microfluidic chamber (diameter <I> = 650/xm, height 
h rsj 80/xm) molded in agarose gel (gray areas) contains the 
lfim janus colloids. H2O2 feeding is achieved by constant cir- 
culation (40/iL/min flow rate) in side channels (dark blue), 
(c) The active colloidal suspension is observed through a 
piezo-driven high numerical aperture objective (Nikon, Wa- 
ter immersion 60x, NA=1.2) mounted on an inverted micro- 
scope. 



crofluidic system based on a gel-microdevice technology 
has been developed, allowing to ensure constant renewal 
of the chemical fuel (H2O2), as well as removal of waste 
products (O2), in an open reactor configuration. On the 
basis of this device, we have explored the behavior of the 
active suspensions at two complementary levels: we have 
first characterized the individual dynamics of active col- 
loids by particle tracking measurements; then we have 
investigated the behavior of a dilute active suspension of 
these particles under an external gravity force field, in 
the same spirit as the historical Jean Perrin experiment 
[2T] . This allowed us to connect the micro-dynamics of 
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individual entities to the macroscopic equilibrium behav- 
ior of the suspension, allowing to probe the link between 
fluctuations and dissipation in this configuration. 

Active colloids and experimental setup - A monolayer 
of commercial fluorescent latex colloids (1 /im diameter, 
Molecular Probes F8823) is formed on a silicon wafer by 
evaporation from 10 3 dilution in isopropanol (99.99%, 
Roth), and coated with 2nm platinum by sputtering, 
see Fig-[l}a. Resuspension is achieved by sonication in 
ultra-pure water (Milli-Q, resistivity 18.2 MOhm.cm -1 ), 
leading, after centrifugation, to 50 /iL of a janus col- 
loidal solution at ~ 0.1% v/v (about 10 9 particles/mL). 
Such janus particles were shown to self-propel in a so- 
lution of hydrogen peroxide, due to the dismutation of 
this chemical on the platinum covering half of the col- 
loids [15,. The properties of this colloidal suspension 
are then investigated in a dedicated microfluidic device, 
sketched in Fig.l-b-c. It is made of a circular microflu- 
idic chamber (diameter $ = 650/im, heights 80/im, vol- 
ume V ~ 30nL) molded in agarose gel [22j[23], and sur- 
rounded by two side channels separated by 125/im gel 
walls. The central chamber is initially filled with janus 
colloids, while H2O2 solution at concentration Co is con- 
tinuously circulated in lateral channels. This gel mi- 
crosystem ensures a constant renewal of H2O2 fuel and 
removal of chemical waste products (O2) by diffusion 
through hydrogel walls from the infinite reservoir and 
sink constitued by the circulating lateral channel. This 
provides a convection-free environnement in the colloids 
chamber with stable chemical conditions over hours, al- 
lowing to study solutions of active colloids (while a sim- 
ilar study would be precluded in a capillary due to the 
production of oxygen bubbles). In this configuration we 
then performed measurements of both the dynamics of 
the colloids and their sedimentation properties, for vari- 
ous concentrations Co of the fuel. 

Colloid micro- dynamics - First, the dynamics of in- 
dividual active colloids is investigated using high speed 
tracking measurements to resolve temporally the differ- 
ent dynamical regimes (see below). For various H2O2 
concentrations Co, the two-dimensional (x,y) motion of 
colloids is recorded with a high speed camera (Phan- 
tom V5) at 100Hz and trajectories are extracted with 
a single particle tracking algorithm (Spot Tracker, Image 
J [24]). The mean square displacement of the colloids 
is obtained as AL 2 (At) = ((R(t + At) - R(t)) 2 ) where 
R(t) is the (2D) instantaneous colloid position and the 
average is performed over time for each individual tra- 
jectory and then over an ensemble of trajectories (typ- 
ically 20). We plot in Fig. [2] AL 2 (t) as a function of 
time for bare (non-active) and active colloids in a so- 
lution of hydrogen peroxide, confirming the impact of 
injected chemical power on the individual motion of the 
colloids, in agreement with Ref. [15]. For the bare (non- 
active) colloids, the dynamics is purely diffusive with a 
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FIG. 2: (a) Experimental mean squared displacements 
AL 2 (A£) and 2D trajectories (inset) for bare (blue) and ac- 
tive colloids (red) in 7.5% H2O2 solution. Bare colloids (bot- 
tom) show standard diffusion (AL 2 linear in time), while 
the mean squared displacement of active colloids is fitted ac- 
cording to Eq.(IT]). The measured diffusion coefficients are 
Do — 0.33/xm 2 /s for bare and D e ff = 1.9/xm 2 /s for active 
colloids, (b) PDF of particle displacement AR for various 
lag time At as a function of the normalized displacement 
AR/AL(At). At=0.3s(*), ls(o), 3s(o) for bare (blue) and ac- 
tive colloids (red) in a solution of 7.5% of H2O2. The dashed 
lines is the gaussian curve. 

diffusion coefficient AL 2 /AAt = D = 0.34 ± 0.02/im 2 /s, 
which is found to be independent of the H2O2 concentra- 
tion. Note that the same value and behavior is obtained 
for janus colloids in the absence of the H2O2 fuel. For 
the janus active colloids in a hydrogen peroxide solution, 
the mean square displacement differs drastically from the 
equilibrium diffusive dynamics and strongly depends on 
the fuel concentration. The colloid exhibits ballistic mo- 
tion at short times, AL 2 (t) ~ V 2 t 2 , while at longer times 
a diffusive regime, AL 2 (t) ~ 4D e fft, is recovered with an 
effective diffusion coefficient D e R much larger than the 
equilibrium coefficient Do- As discussed in pj [15], the 
active colloids are expected to perform a persistent ran- 
dom walk, due to a competition between ballistic motion 
under the locomotive power (with a constant swimming 
velocity V), and angular randomization due to thermal 
rotational Brownian motion. The transition between the 
two regimes occurs at the rotational diffusion time r r of 
the colloids. The characteristic ballistic length scale is 
accordingly a = V x r r . For time scales long compared 
to r r , the active colloids therefore perform a random ran- 
dom walk with an effective diffusion D e R = Do + V 2 r r /4. 
The full expression of the mean squared displacement at 
any time is obtained as [15] : 

AL 2 (At) = ADoAt + ^ [— + e^ - 1] (1) 
2 r r 

We fit the experimental mean squared displacement 
AL 2 (At) using Eq. ([!]) with the propulsion velocity V as 
the only free parameter, while the value of Do is taken 
from the equilibrium diffusion coefficient measured in wa- 
ter and the Stokes expectation is used for r r (r r = 0.9s). 
As shown in Fig. [2}a, an excellent agreement with the 
experimental results is found. Under the present condi- 
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FIG. 3: Normalized density profiles p/po(z) for the active col- 
loidal suspension in stationary state, for increasing swimming 
activity, i.e. increasing D e s. Experimental data (symbols) 
are well fitted by an exponential decay p(z) = po exp(— z/5 e ff ), 
with a sedimentation length £ e ff, which strongly depends on 
(and increases with) the swimming activity. Here £ e ff varies 
from 6pm at equilibrium up to 21/im. 



tions, the measured propulsion velocities V range from 
0.3 pm/s to 2.5pm/ s, corresponding to an increase of 
the effective diffusion coefficient by a factor up to 5.5 
as compared to the bare diffusion. In the following we 
use this measure of D e ^ as a probe of the colloidal 'ac- 
tivity'. Finally we have also measured the probability 
distribution function (PDF) of the colloid displacement 
AR = \R(At) - R(0)) \ for a given time lag At, see Figjij 
b, both for bare (blue) and active (red) colloid particles. 
In both cases, the PDF fits very well to a gaussian with 
variance AL 2 (At) given in Eq. ([!]). Even though depar- 
tures from a gaussian are expected at short time for the 
persistent random walk, these are within experimental 
uncertainty. 

Sedimentation of active particles - We now turn to the 
investigation of the sedimentation of an assembly of such 
active colloids, probed in a sedimentation experiment, in 
the same spirit as the historical Jean Perrin experiment 
[2T] . At thermal equilibrium with a bath at temperature 
T, a dilute population of colloid with (buoyant) mass 
m under gravity g exhibits a steady Boltzmann distri- 
bution profile p{z) = poex.p(—z/So) with So = k^T/mg 
the sedimentation length, which balances gravitational 
and thermal energy. Alternatively this density profile 
can be seen as the stationary solution of the Smolu- 
chowsky diffusion-convection equation, dtp + V • J = 0, 
with J = —D Vp + pmg p the particle flux, and Do? M 
the colloids diffusion coefficient and mobility. This leads 
to the fluctuation-dissipation relationship D = k^Tp. In 
order to explore the validity of these concepts for the 
out-of-equilibrium active suspension, we have measured 
- simultaneously to the individual particle tracking ex- 
periments - the density profiles p(z) of the colloids in 
the microfluidic chamber, for various fuel concentration 
Co. Colloid profiles are measured by scanning the cham- 
ber using a piezo-mounted microscope objective (PIFOC 
P-725.2CD, Physik Instrumente). At each altitude z, 



a stack of images with lateral dimensions 150 x 200/im 
is acquired at 0.3 Hz with a fluorescence camera (Orca, 
Hamamatsu). On each stack, image analysis is performed 
with Mat lab: first the maximum fluorescence intensity 
/foe is determined from "in-focus" colloids, then a 0.3if oc 
threshold criterion is applied in order to discard out-of- 
focus colloids, thus defining a slice with thickness of ± 
2pm [25 . A stack of 100 images is used to obtain a 
good statistical convergence. Overall this allows to ob- 
tain the average number of colloids at each altitude z. 
Note that adsorption of the colloids on the bottom sur- 
face was found to bias the density profiles close to the 
bottom surface and we thus discarded data from the two 
first slices. The results of these " Jean-Perrin" experi- 
ments are presented in Fig. [3] We have checked that a 
stationary state of the sedimentation profile is reached. 
As shown in Fig. |3j the density profiles of the active col- 
loidal suspension p{z) is decreasing with the altitude z 
and, as in the thermal case, can be very well fitted by an 
exponential decay p(z) = po exp(— z/5 e ff)> where po was 
used to normalize the different measurements. The sedi- 
mentation length (5 e ff is found however to depend strongly 
on the activity of the colloids: S e fi increases with an in- 
creased propulsion of the colloids, i.e. injected energy, as 
measured (independently) by their effective diffusion co- 
efficient D e ff. The exponential decay suggests that in the 
present limit of a dilute active suspension, the active col- 
loids still obey an effective Smoluchowsky equation, with 
the current replaced by J = —D e ^Vp + pmg p. This 
predicts a sedimentation length in the form 



^eff = X Ajff? 

Vt 



(2) 



where vt = pmg is the sedimentation velocity in the 
gravity field. 

We have checked this relationship by plotting the sedi- 
mentation length measured from the colloid density pro- 
files, against the effective diffusion coefficient measured 
in the individual tracking measurements. As shown in 
Fig. |1J the predicted proportionnality between sedimen- 
tation length and effective diffuson coefficient is demon- 
strated experimentally for all propelling activities, with 
a proportionality constant which is furthermore found to 
agree with its expected value in Eq. This effec- 

tively connects the micro-dynamics of the active colloids 
to their global stationary profile. 

Finally, we have also measured in the previous experi- 
ments the mean squared number of particles, AN 2 (z) in 
each slice Sz for various altitudes z. We plot in Fig. [4} 
(inset) AiV 2 (z) versus the averaged number N(z) for all 
fuel concentrations. Data are found to collapse onto a 



single curve AN 



• N with 



1. This shows that 



under the present dilute conditions, the active suspension 
behaves like an equilibrium (ideal) solution, and do not 
exhibit anomalously large fluctuations (AN 2 ~ A/" 7 , with 
7 > 1), as previously predicted [1 , however for denser ac- 
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FIG. 4: Sedimentation length £ e fF as a function of effective 
diffusion D e f? extracted from Figs. [2] and [3] respectively. The 
dashed line is a linear fit, £ e ff = a x D e ff as expected from Eq. 
([2]). The measured slope a =19.5 ±2 /im -1 .s is furthermore in 
good agreement with the expected value a — 1/vt — So /Do = 
16.2 zb 2.5 /im -1 .s, with £0 and Do measured independently in 
the absence of injected H2O2. Inset: mean squared number of 
colloids in each slice AN 2 versus the averaged number N for 
various swimming activity. Symbols are similar to Fig [3] The 
dashed line is the prediction for equilibrium system AN 2 = 
k • N, with k = 1 for dilute systems. 



tive systems. 

Our results show that in the present regime, the ac- 
tive colloids behave merely as "hot" colloids, with an 
effective temperature fc^^eS = D^/fi much larger than 
the bare temperature, and do obey an effective Smolu- 
chowsky equation. In the Wun and tumble^ model of 
particles under external fields discussed by Tailleur and 
Cates, this behavior is expected in the regime where the 
swim speed V is larger than the sedimentation velocity 
vt- Here V ~ /im/s, while vt = 1/ck ~ 50.10 _3 /im/s 
- see Fig. ^ , and one has V/vt ^> 1 in the present 
experiments. Our experimental results thus validate the 
theoretical expectations in this regime. 

To conclude, we have proposed efficient experimental 
tools allowing to explore the properties of active col- 
loidal suspensions under controlled and tunable condi- 
tions. This involves artificial chemically active colloids, 
studied in a dedicated gel microffuidic system. This 
opens up many perspectives towards a thorough explo- 
ration of the out-of-equilibrium behavior of active supen- 
sions, in particular for higher volume fractions, where 
exotic behavior was predicted theoretically, but not ob- 
served up to now with artificial active fluids. This in- 
volves collective effects, the emergence of flow order, 
anomalous fluctuations [J, nematic ordering [7], etc. But 
beyond these predicted behaviors, it is interesting to note 
that the present colloids not only interact via hydrody- 
namic flows, but also via chemical interactions, through 



the spatial (diffusive) extension of the consummed fuel. 
This is expected to affect the many-body behavior of the 
suspension, with couplings which have not been consid- 
ered up to now in the litterature. Work along these lines 
is in progress. 
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